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^ . We study the two point correlation function of two nucleon currents 

sandwiched between the vacuum and the pion state. The light cone 
QCD sum rules are derived for the ttNN coupling g n NN an d the 
7riV7V(1535) coupling g n NN*- The contribution from the excited 
states and the continuum is subtracted cleanly through the double 
£T) ', Borel transform with respect to the two external momenta, p\-,p\ = 

(p — q) 2 . We first improve the original sum rule for g-^NN and 
determine the value of the pion wave function ^^{uq — |) = 1.5 ± 
0.2, which is a fundamental nonperturbative quantity like the quark 
condensate. Our calculation shows that the 7riViV(1535) coupling 
is strongly suppressed. 

PACS Indices: 13.75.Gx; 14.20.Gk; 14.40.Aq; 13.75. Cs; 12.38.Lg 
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1 Introduction 



The nNN coupling g^NN and the 7riViV(1535) coupling g-^NN* play a very important role 
in one boson exchange potentials (OBEP) for the nuclear force. Although it is widely 
accepted that QCD is the underlying theory of the strong interaction, the self-interaction 
of the gluons causes the infrared behavior and the vacuum of QCD highly nontrivial. 
In the typical hadronic scale QCD is nonperturbative which makes the first principle 
calculation g n NN and g n NN* unrealistic except the lattice QCD approach, which is very 
computer time consuming. So a quantitative calculation of the nNN and 7riViV(1535) 
with a tractable and reliable theoretical approach proves valuable. 

The method of QCD sum rules (QSR), as proposed originally by Shifman, Vainshtein, 
and Zakharov [|TJ and adopted, or extended, by many others ^ f|, are very useful in 
extracting the low-lying hadron masses and couplings. In the QCD sum rule approach the 
nonperturbative QCD effects are partly taken into account through various condensates 
in the nontrivial QCD vacuum. In this work we shall use the light cone QCD sum rules 
(LCQSR) to calculate the nNN and vrAW(1535) couplings. 

The LCQSR is quite different from the conventional QSR, which is based on the short- 
distance operator product expansion. The LCQSR is based on the OPE on the light cone, 
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which is the expansion over the twists of the operators. The main contribution comes from 
the lowest twist operator. Matrix elements of nonlocal operators sandwiched between a 
hadronic state and the vacuum defines the hadron wave functions. When the LCQSR 
is used to calculate the coupling constant, the double Borel transformation is always 
invoked so that the excited states and the continuum contribution can be treated quite 
nicely. Moreover, the final sum rule depends only on the value of the hadron wave function 
at a specific point, which is much better known than the whole wave function 0]. In the 
present case our sum rules involve with the pion wave function (PWF) (p n (uo = \)- Note 
this parameter is universal in all processes at a given scale. In this respect, (p n (uo = |) is 
a fundamental quantity like the quark condensate, which is to be determined with various 
nonperturbative methods. Like the quark condensate, it can be determined through the 
analysis of the light cone sum rules. In |J the value is obtained as (p n (uo = \) = 1.2 ±0.3 
using the pion nucleon coupling constant and the phenomenological pum coupling constant 
as inputs. 

The LCQSR has been widely used to derive the couplings of pions with heavy mesons 
in full QCD ||, in the limit of rriq — > oo @ and 1/tuq correction [jlO[| , the couplings 



of pions with heavy baryons ||11|| , the p decay widths of excited heavy mesons fi~2"f , and 
various semileptonic decays of heavy mesons |L3| etc. 

The aim of the paper is two- fold. First we improve the original sum rule for g^NNiq 2 — 
0) || in the following way: (1) the contribution from the gluon condensate < g 2 s G 2 > 
and the quark gluon mixed condensate (g c qcr ■ Gq) are calculated, which is numerically 
not negligible as in the case of the nucleon mass sum rule || and the sum rule for the 
pion nucleon coupling constant |I3fl ; (2) the twist four contribution is collected in a more 
transparent form and is estimated slightly differently from ||; (3) the uncertainty due 
to Ajv is reduced in the numerical analysis with the help of the Ioffe's mass sum rule. 
We arrive at (p n (l/2) = 1.5 ± 0.2 using the experimentally precisely known g n NN |T5 ] . 



Secondly we employ the LCQSR to calculate 7riViV(1535) coupling. The continuum and 
excited states contribution is subtracted rather cleanly within our approach. Our result 
shows that the 7riViV(1535) coupling is strongly suppressed. In both cases the final sum 
rules are stable with reasonable variations of the Borel parameter and the continuum 
threshold. 

Our paper is organized as follows: Section [l] is an introduction. We introduce the two 
point function for the 7vNN vertex and saturate it with nucleon intermediate states in 
section 0. The definitions of the PWFs are presented in section |3|. In the following section 
we present the LCQSR for the ttNN coupling and discussions of these PWFs. In section 
|5] we employ LCQSR to calculate the 7riViV(1535) coupling. A short summary is given in 
the last section. 



2 Two Point Correlation Function for the 7tNN coupling 

The pion nucleon coupling constant g n NN has been calculated with the following variations 
of the traditional QCD sum rule method: (1) the three point correlator of two nucleon 
and one pseudo-scalar meson interpolating fields, which are saturated with resonances in 
the nucleon N and pion tt channels on the phenomenological side 0, [T7[ |TBJ; (2) the 
two point function of two nucleon interpolating fields sandwiched between the vacuum 
and one 7r state and saturating only with N resonances 0, [T7|, [19]; (3) the external field 
method, which considers the two point correlator of two nucleon interpolating fields in 
the presence of the pion field |T^| ; (4) the light cone QCD sum rules ||; (5) using the 
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Goldberger-Treiman relation with the nucleoli axial coupling constant determined by 
QCD sum rules in external fields [0. The recent calculation of g& = 1.37 ± 0.10 H yields 

4™ = 14.7 ±2.1. 

4-7T 

In the method (1) the singular pole structure ^pr was picked out and identified for 
the extraction of g n NN- Note the operator product expansion (OPE) for the correlator 
(0) is valid only in the region p\ <C 0, p\ < 0, q 2 < 0. At q 2 = the OPE does not 
hold. Moreover it was pointed out that with the first methods the excited states and 
continuum contribution was not subtracted away |3], 20], which could contaminate the 
sum rule severely. 

We start with the two point function 



n(pi,p 2 ,9) = J d 4 xe^ x {0\T Vp (x)ff n {0)\^+(q) i 

with pi = p, p 2 = p — q and the Ioffe nucleon interpolating field || 

n P {x) = e abc [u a (x)C lfl u b (x)] l5 Yd c (x) , (2) 

fjpiv) = e abc [u b (yh u Cu aT (y)}d c (y)Yl 5 , (3) 
where a, b, c is the color indices and C = ^7270 is the charge conjugation matrix. For the 
neutron interpolating field, u «-> d. 
n(pi, P2, q) has the general form 

n(pi,p 2 , q) = F(pl,Pl, Q 2 )Ql5 + Fi(p 2 ,pi, q 2 h 5 + F 2 (p 2 v pi, q 2 )Pl 5 + F 3 (pl, p 2 2 , q 2 )o llv ^p il q v 

(4) 

It was well known that the sum rules derived from the chiral odd tensor structure yield 
better results than those from the chiral even ones in the QSR analysis of the nucleon 
mass and magnetic moment @, ^lj. Most of the QSR analysis of the pion nucleon coupling 
constant deals with the tensor structure 575. Based on these considerations we shall focus 
on the same tensor structure and study the function F^p^p 2 ., q 2 ) in detail. 

The pion nucleon coupling constant g n NN is defined by the 7riV interaction: 

ZttNN = g^NNNi^r ■ ttN + h.c. , (5) 

where h.c. stands for the hermitian conjugate. 

At the phenomenological level the eq.(|36|) can be expressed as: 



n( Pl)P2) q) = ^m^^(g 2 ) (pf _ M ^_ M ^ + • • • (6) 

where we include only the tensor structure 75^ only. The ellipse denotes the continuum 
and the single pole excited states to nucleon transition contribution. A at is the overlapping 
amplitude of the interpolating current T)n(x) with the nucleon state 

(0|77(0)|JV(p)) = A^(p) (7) 

3 The Formalism of LCQSR and Pion Wave Functions 

Neglecting the four particle component of the pion wave function, the expression for 
F(p 2 ,p2,q 2 ) with the tensor structure at the quark level reads, 

J e^dx(0\T Vp (x)n n m^(q)) = ^e abc e a ' b ' c '^ 

(8) 
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where iS(x) is the full light quark propagator with both perturbative term and contribu- 
tion from vacuum fields. 



iS(x) = (0\T[q(x),q{0)]\0) 
x (qq) x 2 



2n 2 x 4 12 192 



(qg s a ■ Gq) (9) 



lD7T^ JO X 1 X 1 

where we have introduced x = x^ . The relevant feynman diagrams are presented in 
FIG 1. The squares denote the pion wave function (PWF). The broken solid line, broken 
curly line and a broken solid line with a curly line attached in the middle stands for the 
quark condensate, gluon condensate and quark gluon mixed condensate respectively. 

By the operator expansion on the light-cone the matrix element of the nonlocal op- 
erators between the vacuum and pion state defines the two particle pion wave function. 
Up to twist four the Dirac components of this wave function can be written as ||: 

i 

lu 1 x (,n (o,\ J- rr^n- (o,\ _L n( 



<7r(g)|rf(x) 7/ , 75 n(0)|0> = -if„qj du e luqx (<pM + x 2 9l (u) + 0(x 4 )) 

jo 

+ U{x, - ^) / du e iuqx g 2 (u) , (10) 

QX JO 



< n(q)\d(x)i-yr oU (0)\0 > = C du e iuqx ip P (u) , (11) 



m u + m d Jo 

— f Tfl 2 

< n{q)\d{x)a^ 5 u{0)\0 > = i{q^x u - q u x^) Jw n / du e luqx ^{u) . (12) 



< ■K(q)\d(x)a a pj 5 g s G IJiU (ux)u(0)\0 >= 



and 



< n(q)\d(x)-f^ 5 g s G a(3 (vx)u(0)\0 >= 

xpq» 

q ■ x 



U[qp{pan - -^ Ji ) - q a (. 



q ■ x 



90n 



Vaitp ± (ai)e 



iqx(a.\-\-va.z) 



Q 

+fn—^-(q a Xi3 - qpx a ) I Va.iip\\(ai)e 



q ■ x 



iqx(cti+va 3 ) 



(14) 



< Tr(q)\d(x)y lJl g a G a p(vx)u(0)\0 > 



9afi 



Xa.q\i 

q ■ x 



q a \9fiii 



xpq^ 

q ■ x 



iqx(ai+vas) 



q ■ x 



(q a xp - qpx a ) I Vaiip\\(a.i)e 



D iqx{a\+votz) 



(15) 



The operator G a p is the dual of G a p\ G a /3 = ^ai3S P G 5p ; Vcti is defined as Voti = 
da\da2da^5{l — a± — a 2 — ot-s) ■ Due to the choice of the gauge x^A^{x) = 0, the path- 
ordered gauge factor Pexp (ig s Jq 1 dux^ A^ux)) has been omitted. The coefficient in front 
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of the r.h.s. of eqs. ([H]), (|T2"D can be written in terms of the light quark condensate 

m 2 2 

< uu > using the PCAC relation: fx n = — = — — < uu >. 

m u + m d f£ 

The PWF ip n (u) is associated with the leading twist two operator, gi(u) and g 2 (u) 
correspond to twist four operators, and y>p(u) and (p a (u) to twist three ones. The function 
y? 37r is of twist three, while all the PWFs appearing in eqs. ([13]), ( |T5D are of twist four. The 
PWFs ip(xi, fx) (fx is the renormalization point) describe the distribution in longitudinal 
momenta inside the pion, the parameters x i — 1) representing the fractions of the 

longitudinal momentum carried by the quark, the antiquark and gluon. 

The wave function normalizations immediately follow from the definitions (fL0|) -(|l5D: 
f^du (Ptt(u) = Jodu ip a (u) = 1, Jodu gi(u) = S 2 /12, fVaiip±(ai) = f Vaii(p\\(a£i) = 
0, / Ucti(p±(ai) = — J Uati<p\\(ai) = 6 2 /3, with the parameter 6 defined by the matrix 
element: < n(q)\dg s G a ^'y a u\0 >= i5 2 f^q^. 



4 The LCQSR for the ttNN coupling 
Expressing @ with the pion wave functions, we arrive at: 

U( Pl ,p 2 ,q) = I d'x [ 1 due i ^>{^{[^(u)+x 2 g l (u)} l5 q-i l5 (x-— qg 2 (u)} 
--^H^{qq)Va(u)l b (x 2 q -q-xx)- Jl fx n ml(qq) ^ a (u)^(x 2 q - q ■ xx) 

-i J d A x £du J Pa 4 e^-^ Ql+u ^^ 

+2(1 - 2u)(q ■ x)q] + ■ ■ ■ ,(16) 

where fx n = 1.65GeV, f m = 132MeV, (qq) = -(225MeV) 3 , (qg s a ■ Gq) = m 2 (qq), m 2 = 
0.8GeV 2 , a = — (2Tr) 2 (qq). We have collected the terms relevant for the tensor structure 
7 5 g in ([TJ) only. 

In the following sections we will frequently use integration by parts to absorb the 
factors (q ■ x) and l/(q ■ x), which leads to the derivative and integration of PWFs. For 
example, 



q ■ x 



^(u)e-^ x du = -i —-^du + ^(u)e-™«-% , (17) 



o (x 2 ) n 

„l p -iuq-x „ 

/ g 2 (u)du = -i / e~ iuq - x G 2 (u)du - G 2 (u)e 

Jo q ■ x J 

where the functions G 2 (u) is defined as: 



iuqx 1 1 
10 ' 



G 2 (u) = - / g 2 (u)du . (19) 
Jo 

Note the second term in (|I7|) and fll8D vanishes due to (p n (uo) = G 2 (uq) = at end points 
m = 0, 1. 

We first finish Fourier transformation in ([16]). The formulas are: 



f ^__ d n x = z( _ ir+1 * D - 2n - D/2 r(D/2-n) 

J (X 2 Y H ' f-p2)D/2-n T(n) > ^> 
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*e*> T _ f ir+1 2 g - 2 " +1 vr^ 2 r(D/2 + l-n) 

(x 2 )« 1 J (_ p 2 ) D/2+l-n r(n) P- W 

Making double Borel transformation with the variables p\ and p\ the single-pole terms 
in (|6|) are eliminated. The formula reads: 

^I ^fI^ = (MV-»e^^ - «„) ■ (22) 

Subtracting the continuum contribution which is modeled by the dispersion integral 
in the region Si, s 2 > So, we arrive at: 

V2m N X N g 7TNN e ~m* = 
■{+ A^ K)M 6 /2( _^) _ ^[ gi ( Uo ) + G 2 (n )]M 4 /i(^) 



7"G/„, \ i tG(„, \ otG(„, m i/Stt 2„, rG/ 



« M /o(th)^W + i 2 K) - 2J 3 G ( Mo )] - ^W« Q i 4 >o)} , (23) 



6vr 2 JUV M 2n 1 v u/ 2 v u/ 6 v U/J 6tt 



it 



where f n (x) = 1 — e 11 J] fr- is the factor used to subtract the continuum, sq is the 

fc=0 

continuum threshold. «o = M -f+ M 'i , ^i 2 = m'^+m' 1 ' ^i ' ^2 are the Borel parameters, and 
The functions lf(ua) are defined as: 

if (no) = f ^ 1 ^ (ai ' 1 ~" '"°" ai) , (24) 

JO Mo — Q!l 

i 2 (tt ) = / cto 2 : , (25) 

JO 1 — «o — ^2 

Jfta) = f da, da /3,( ai ,a 2 ,l- ai - a2 ) (26) 

JO JO (1 — Qfl — a 2 r 



r G/ \ p. f^" ^ ^ «i,a 2 ,l -ai -a 2 

if ( u o) = / aai / da 2 (1 - 2w + «i - a 2 )«o ■ (27) 

Jo jo (1 — oli — a 2 ) 

The twist four PWFs ip±(aii), <P\\( a i), l P±( a i) an d <P\\(ck) do not contribute to the 
chiral odd tensor structures, which was first observed in 0. Moreover our sum rule ( |23| ) 
is symmetric with the Borel parameters M 2 and Mf . So it's natural to adopt M\ = Mf = 
2M 2 , m = 2 - We shall work in the physical limit q 2 = m 2 — > and discard the terms 
with the factor q 2 in (p3|) . 

The various parameters which we adopt are a = 0.55 GeV 3 , (g 2 G 2 ) = 0.48 GeV 4 , 
m 2 = 0.8 GeV 2 at the scale /i = lGeV, s = 2.25GeV 2 , m N = 0.938GeV, X N = 0.026GeV 3 

!• i 

The working interval for analyzing the QCD sum rule (H) is 0.9GeV 2 < M| < 
1.8GeV 2 , a standard choice for analyzing the various QCD sum rules associated with the 



G 



nucleoli. In order to diminish the uncertainty due to A at, we shall divide (|23|) by the 
famous Ioffe's mass sum rule for the nucleon: 

32^-3 = *•/,<£) + \ M ^) + \«> ~ ^ • (28, 

The resulting sum rule depends on the PWFs, the integrals and derivatives of them 
at the point uo = |. Since 5 2 is numerically small, the uncertainty due to the integral 
term G^wo) is insignificant. There are many discussions about the leading twist PWF 
(p n (u) in literature [2^, ^3], |, [24], 25], [26], 27J. For example, at uq — \ the values of various 



model functions are: ip n (u ) = 1.22 §, 1.273 1.25 [g|], 1.71 |24|, 1.35 |27| and 1.5 
for the asymptotic PWF respectively. 

Knowledge of the PWFs involved with the gluon field is still very limited. Only the 
very lowest a few moments of these PWFs were calculated with significant errors using the 
method of QSR, which in turn were used to determine the detailed forms of the functions. 
Such an approach is sometimes misleading since there are many wave functions satisfying 
the same constraints from moments. 

In order to illustrate this point more clearly we use the determination of <p n (u) as an 
example. The model wave function for <p n (u) based on the QCD sum rule approach was 
given in [[J as: 

/ 3 15 

V?7t(m, fi) = Quu ( 1 + a2(/u)-[5(w — u) 2 — 1] + a 4 (/i) — [21(w — u) A — 14(w — u) 2 + 1 
V 2 8 

(29) 

where u = 1 — u. Yet the authors pointed out themselves that the oscillations of (|29|) 
around uq = | is unphysical, which is due to the truncation of the series and keeping only 



the first a few terms when ^(w) is expanded over Gegenbauer polynomials. In [25] it 
was stressed that: (1) the expansion over Gegenbauer polynomials converges very slowly 
as can be seen from the large value of a 2 and a^, (2) any oscillating wave function is not 
physical, since no detached scale is seen to govern such oscillations. Recently Mikhailov 
and Radyushkin reanalyzed the QCD sum rules for ^(m) taking into account the 
non-locality of the condensates. They suggested the following wave function: 

g 

<Ptv{u) = —Vuti . (30) 

7T 

Halperin suggested the following form [p5] : 

(p v (u) = Nexp I --^r- ) [fi 2 +w+ - 2)uu) , (31) 




where N = 4.53, n = ||, jl = ^ Urn* + fm p J = || with m = 330MeV, m = 620MeV, 
and j3 = 320MeV. Note all the above three forms of PWFs are rather close to the asymp- 
totic form 

^asym (u) = gu - ? (32) 

which are in strong contrast with the original double humped Chernyak-Zhitnitsky form 
keeping the lowest two orders of the expansion f22[: 

y^ sym (w) = 30mm(1 - Auu) . (33) 

In other words, after summing the whole series the physical PWF can not deviate too 
much from the asymptotic form. Based on the above consideration we use the asymptotic 



forms for the PWFs involved with gluons as in |TT . 
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For the other PWFs we use the results given in || since they are relatively better 
known. ip P (u ) = 1.19, (p a (u ) = 1.47, gi(u ) = 0.022GeV 2 , ^2(^0) = and Gi{uq) = 
0.02GeV 2 at u — ~ and /i = lGeV. Their first derivatives satisfy: v^K^o) — V^pl^o) = 
^(«o) = g'i(u ) = G' 2 (u ) = 0. 

The dependence on the Borel parameter M 2 of g n NN are shown in FIG 2 with so = 2.25 
GeV 2 using different values of (p n (uo). From top to bottom the curves corresponds to 
y?7r(^o) — 1-6, 1.5, 1.4 respectively. The final sum rule is very stable in the working region 
of the Borel parameter M 2 and sensitive to the value of (p n (^)- Experimentally the tcNN 
coupling constant has been extracted very precisely: g^NN = 13.5 [|l^]. Using this value 
as the input we obtain: 

^(m ,m= lGeV) = 1.5 ±0.2 , (34) 

which is very close to the asymptotic PWF. 

It is interesting to notice that the QSR derived from the chiral odd structure 975 
yielded (p n (u Q ) = 1.6 while the QSR from the chiral even structure qfry 5 led to a smaller 
value 0. After averaging the two values the authors got (p n (u ) = 1.2 ± 0.3. As pointed 
out in Section |2], generally the QSR derived from chiral odd structure is more reliable. 
Therefore, such an averaging may not be very suitable here. 

Very recently Belyaev and Johnson investigated the relation between light-front quark 
model and QCD |27j. They found additional support that the two point PWF ^^{u) at- 



tains the asymptotic form. Numerically <p n (u) at u — | agrees with the asymptotic PWF 
within 10%, which was consistent with our result. It appears that almost all the phe- 
nomenological analyses share the common feature, i.e., that the PWF starts to approach 
the asymptotic form more or less at the low scale [i = lGeV already. This point deserves 
further investigation. 



5 The LCQSR for the ttAW(1535) coupling 

Quantum Chromodynamics (QCD) is asymptotically free and its high energy behavior has 
been tested to one-loop accuracy. On the other hand, the low-energy behavior has become 
a very active research field in the past years. Various hadronic resonances act as suitable 
labs for exploring the nonperturbative QCD. The inner structure of nucleon and mesons 
and their interactions is of central importance in nuclear and particle physics. Interna- 
tionally there are a number of experimental collaborations, like TJNAL (former CEBAF), 
COSY, ELSA (Bonn), MAMI (Mainz) and Spring8 (Japan), which will extensively study 
the excitation of higher nucleon resonances. 

Among various baryon resonances, negative parity resonance A^*(1535) is particularly 
interesting, which dominates the 77 meson photo- or electro-production on a nucleon. 
The branching ratio for the decay N* — > r]N is comparable with that for N* — > ttN. 
Considering the phase space difference and using the experimental decay width of A"* 
p8| , we get g^N* = 2 and g n NN* = 0.7. The latter value is in strong contrast with the 
pion nucleon coupling g n NN = 13.4. Thus arises the question: why is the coupling g n NN* 
so small compared with g n NN^ 

Whether the coupling g n NN* is strongly suppressed is under heated debate in litera- 
ture. In a recent extended coupled channel analysis of irN scattering, the Jiilich group 
used g V NN* = 1-94 and g^NN* < 0.12 Jido, Oka and Hosaka suggested in a recent 
letter that ttNN* coupling is strongly suppressed as a consequence of chiral symmetry 
3T!]. Their argument was based on a pion-nucleon correlator of two baryon interpolating 



fields. The chiral transformation properties of their interpolating fields then imply that 
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the correlator is purely proportional to the tensor structure 75, with no piece of the form 
P7 5 , which is the relevant structure for ttNN* coupling. Based on the above observation 
they claimed that the coupling ttNN* vanishes. The above argument was criticized by 



Birse [31]. He pointed out that the ^75 piece of the correlator is a sum of all possible pion- 
baryon couplings that can contribute. Hence the absence of a P75 piece is a statement 
about the particular combination of the pion-baryon coupling and the subtraction terms 
that correspond to the chosen interpolating fields. It does not imply that the physical 
ttNN* coupling is suppressed. With an interpolating field with covariant derivative for 



the N* Kim and Lee used QCD sum rules to estimate cj-knn* ~ 1.5 f32|. But in their anal- 



ysis the continuum and excited states contribution is poorly subtracted and the numerical 
results depend strongly on the value of the quark gluon mix condensate (qg s o ■ Gq), which 
renders their conclusion less convincing. 

In this section we shall employ the LCQSR to calculate TtNN* coupling. The con- 
tinuum and excited states contribution is subtracted rather cleanly within our approach 
through double Borel transformation. 

We start with the two point function 

n(pi,p 2 ,<?) = / d^xe^ (0\T Vp (x;s)fr n (0;t)\7r + (q)) (35) 

with pi = p, P2 = p — q and the general nucleon interpolating field without derivatives 
which couples to both positive and negative parity nucleon resonances 

r] p (x; s) = e abc { [u a {x)Cd\x)\ j 5 u c {x) + s [u a {x)C^d b {x)\ u c {x)} (36) 

where s is the mixing parameter and rf p = Tjp'jo. Note we use a slightly different nucleon 
interpolating field ( |3"6| ) in this section. The current r) p (x; s = —1) is one half of the one 
defined in Q. Both of them are called Ioffe's current. The later form appears more in 
literature. The Ioffe's current couples strongly to the positive parity nucleon || while 
it was found that the current rj p (x; s = 0.8) is optimized for negative parity nucleons and 
couples strongly to iV(1535) P3J. 

It is important to note that the diagonal transitions like N — > N, N* — ► iV* does 
not contribute to the tensor structure ^75. In other words, the function F 2 involves solely 
with the process iV* — > N and corresponding continuum contribution. Based on the above 
observation we shall focus on the chiral odd tensor structure pjs in this section. 

The 7riV./V(1535) coupling constant g n NN* is defined as: 

ZttNN* = 9ttNN*Nt ■ tcN + h.c. (37) 
At the phenomenological level the eq.(|35|) can be expressed as: 



n(pi,p 2 ,g) = -(mjv+mtfOftrJVJviT-a ^TTT-— ^2 \( 2 iT?2T/ z P75+- 



X N (s)X N *(t) \ N *(s)\ N (t) 

{pl-Ml)(pl-Ml*) (pl-M 2 N *M-M* NJ 

(38) 

where we write the structure py*, explicitly only and the continuum contribution is denoted 
by the ellipse. Xn(s) is the overlapping amplitude of the interpolating current t]n(x) with 
the nucleon state 

(0|77(0; s)\N(p)) = X N (s)u N (p) (39) 

We neglect the four particle component of the PWF and express (p5|) with the PWFs. 
After Fourier transformation and making double Borel transformation with the variables 



9 



p 2 and p\ the single-pole terms are eliminated and finally we arrive at: 



(mjy + mAr»)[AAf(s)AAf*(t) - \N*(s)\ N (t)]g nNN *e M i M 2 



JV ;[7(1 + st) + 3(s + «)] 9 2 («o)M 4 /i(^) + ^(1 + ^ + * - 38t)a^^(« )M 2 /o( " 1 



327r 2L-V "V "V" VJ^V-uy ^V^/ 2 g 87r2 V - / WctV -u, ^Vj^ 

f 46oi^ (1 + S){1 + ^sG^'A^M 2 /^) - ^^(5 + 3s + 3* - llst)am 2 ^M 
F 19^ (S " *)amjAWp(«o) - [19(1 + sf) + 7(a + t)](^G 2 )^ 2 ( Mo ) 



/t 'l + s )(l + i)(s ( J GViW (40) 



11527T 2 



The twist three PWF <^ 37r appears in the combination /Wy^ (<?<?) <? 2 - In the physical 
limit q 2 = m 2 —>■ 0, its contribution is negligible, which is in contrast with the sum rule 
for g^NN (p3|), where the twist three PWF (p 3n plays an important role. 

It is a common practice to adopt the Ioffe current for the ground-state nucleon ||, 
i.e., letting s = —1. Now the equation (pEOP has a simple form: 

(m^v + m N *)[\ N (s)\ N *(t) - A Af »(s)A A r(t)]5f 7rA r 7 v*e M i M 2 = 



k 2V ^ v y J1V M 2 ' 72vr 2 uy JUV M 2 

rl^Q + 7t)aml^ip' a (uQ) - in {^ 2 (l + t)amlfj^(pp(u ) 



3456vr 2V 7 192^ 

A :i-t)^ 2 G 2 )^K) , (4i) 



1927r' 
where t = 0.8. 

The various parameters are si = 3.2GeV 2 , m^* = 1.535GeV, Atv(s = —1) = 0.013GeV 3 
X N *(s = 0.8) = 0.27GeV 3 with the formulas in 0. Moreover X N (s = -1) > X N (t = 0.8), 
\ N *(t = 0.8) ^> \n* (s = —1), so it is reasonable to discard the term Atv(£ = 0.8)Aat*(s = 
— 1) in the eq. (O). We use the model PWFs presented in || to make the numerical 
analysis. 

Our sum rule (|4T| ) is asymmetric with the Borel parameters Mf and Mf due to 
the significant mass difference of N and iV(1535). It is natural to let Mf = 2m 2 N f3, 
Mf = 2m 2 N *f3, where /3 is a scale factor ranging from 1.0 to 2.0. In this case we have 
M 2 = 1.28/3GeV 2 , u = 0.27, g 2 (u ) = -0.03GeV 2 , ^ P {u Q ) = 0.66 and ^(u ) = 2.63 at 
the scale \i = lGeV. 

The sum rule for g^NN* is very stable with reasonable variations of s\ and M 2 as can 
be seen in FIG 3. Finally we have 

9kNN* = {-) (0.08 ±0.03). (42) 



Our result is consistent with the conclusions in Refs. |3C 

We have included the uncertainty due to the variation of the continuum threshold 
and the Borel parameter (3 in fl42f) . In other words, only the errors arising from numerical 
analysis of the sum rule flllD are considered. Other sources of uncertainty include: (1) 
the truncation of OPE on the light cone and keeping only the lowest twist operators. 
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For example the four particle component of PWF is discarded explicitly; (2) the inherent 
uncertainty due to the detailed shape of PWFs in different models; (3) throwing away 
the term Xwit = 0.8)\^*(s = —1) in the eq. (|4T|); (4) the continuum model etc. In the 
present case the major error comes from the uncertainty of PWFs since our final sum 



rule (|¥TD depends both on the value of PWFs and their derivatives at u$ = 0.27. Luckily 
around uq = 0.27 different model PWFs have a shape roughly consistent with each other. 
So a more conservative estimate is to enlarge the error by a factor of two. Now we have 

9*nn* = (-)(0.08±0.06) . (43) 

We want to point out that one should not be too serious about the specific number. 
What's important is the fact that \g n NN* \ |fl , 7r jvjv|- 



6 Discussion 



In summary we have constrained the parameter using the experimentally precisely 

known g n NN using LCQSR. We also have also calculated the tiNN* coupling. The contin- 
uum and the excited states contribution is subtracted rather cleanly through the double 
Borel transformation in both cases. Our result shows explicitly the suppression of g n NN*- 
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ure Captions 



FIG 1. The relevant feynman diagrams for the derivation of the LCQSR for nNN and 
7r7ViV(1535) coupling. The squares denote the pion wave function (PWF). The broken solid 
line, broken curly line and a broken solid line with a curly line attached in the middle stands 
for the quark condensate, gluon condensate and quark gluon mixed condensate respectively. 

FIG 2. The sum rule for g^NN as a functions of the Borel parameter M 2 for different <fiir(^) 
with the continuum threshold so = 2.25GeV 2 . From bottom to top the curves correspond 
to <£7r(§) = 1.6, 1.5, 1.4 respectively. 

FIG 3. The sum rule for g^NN* as a function of the scale parameter (3 with the continuum 
threshold si = 3.4, 3.2, 3.0 GeV 2 using the PWFs in |]]. 
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